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We extend the circle theorem on the zeros of the partition function to a 

^" 

0^ ■ continuum system. We also calculate the exact zeros of the partition function 



for a finite system where the probability distribution for the order parameter 



i-rt ■ is given by two asymmetric Gaussian peaks. For the temperature driven first 

a: 

order transition in the thermodynamic limit, the locus and the angular density 



of zeros are given by r = e^ Ac ' 2l ' e and 2irg(9) = l(l + ^(Ac/l) 2 9 2 ) respectively 

in the complex z(= re i6l )-plane where I is the reduced latent heat, Ac is the 

discontinuity in the reduced specific heat and z = exp(l — T c /T). 
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One of the fascinating subjects of equilibrium statistical mechanics is to understand how 
an analytic partition function acquires a singularity when the system undergoes a phase 
transition. For the last three decades main focus has been on the second order transition. 
Only recently a renewed interest in the first order phase transition began to emerge |L|. 

Since Yang and Lee first published their celebrated papers on the theory of phase 
transitions and the circle theorem on the zeros of the partition function, there have been 
many attempts to generalize the theorem ||. Fisher || initiated the study of zeros of the 
partition function in the complex temperature plane, and Jones f| proposed a scenario for 
the first order transition for a continuum system. However very little is known about the 
distribution of zeros for the continuum case. This is because the partition function for the 
continuum system is not a polynomial in general and the original proof of the circle theorem 
relied heavily on particular properties of the coefficients of a polynomial. Recently we have 
been able to prove the theorem in a quite different approach || and this approach allows us 
to extend the theorem to the continuum case. 

We found that the circle theorem follows from a certain mathematical relation which 
exists between a probability density function and the zeros of its characteristic function. In 
this paper we first prove three theorems and a corollary without referring to the partition 
function. When these theorems and the corollary are translated in physical terms we find 
that, (1) the zeros of the partition function can be expressed in terms of the discontinuities 
in the derivatives of the free energy across the phase boundary if there is a non-vanishing 
discontinuity in the first order derivative; (2) there is no zeros in the single phase region where 
the probability distribution is given by a single Gaussian peak; (3) the zeros of the partition 
function are calculable exactly at the two phase coexistence point where the probability 
distribution is given by two asymmetric Gaussian peaks; and (4) the zeros lie on the unit 
circle if the transition is symmetric. 

Furthermore we find the finite size scaling very much similar to that of the discrete 
system 0. Therefore this result can again be used, for a continuum system, (1) to resolve 
the recent controversy over equal weight vs. equal height of the probability distribution 



functions |7|-f|; and (2) distinguish the first order transition from the second fllOlfl, just as 
we have done for the discrete system in ref. |J. 

Consider a probability density function f(x) of a random variable X of continuous type, 



which satisfy i) f(x) > and ii) J^° 00 f(x)dx = 1. The characteristic function [[D]] of a 
random variable X is defined by 

/oo 
e™f{x)dx (1) 

-co 

where i is the imaginary unit. The logarithm of (j)(u) is known as the second characteristic 
function or cumulant generating function and denoted by ip{uj). That is, 

00 (iuoY 
#,) = ln(0H) = £ 7s ^U (2) 

8=1 S - 

The expansion coefficients, 7 S are known s-th cumulants or semi-invariants and are calculable 
by the formula, 7 S = d s ip(u})/d(iu) s \ iu)=0 . 

If 4>(u) can be analytically continued into the complex cu-plane, M(t) = <p(—it) is a real 
function called the moment generating function. Since the characteristic function always 



exists and its properties are well known [IT|, we will consider zeros of the characteristic 
function. 

Let S* be a subset of the sample space of a random variable X. The conditional distri- 
bution of X within the sample space, S 1 *, is called the truncated distribution of X, whose 
probability density function /*(x) is given by 

M ' )= i$k- (3) 

where ie5,. Suppose the characteristic function of the distribution which can be analyti- 
cally continued into the complex plane. We will consider zeros of the characteristic function 
in the finite region only. If the density function is not a delta function (in that case there are 
no zeros in the finite region) we can always divide the density function into two truncated 
density functions, taking a point x* some point in the middle of the distribution. 
That is, 



f( x ) = c (/i(x) +af 2 {x)), 



(4) 



where 



/i(*) 



/(*) 



JTL /w^ 



>Aj "\ >Aj ;}; 



: a; > a;* , 



/*(*) 







/W 



c = I-oo f(x)dx, and a = f™ f(x)dx/ /f ^ /(x)dx. 
Consider the zeros of characteristic function 

0(a;) = 0i (a;) + a</> 2 (a;). 



*/y -^ tC' -fe 



(5) 



If V'iX 1 - ') an d ^2 (<*>), the cumulant generating functions, of 4>i(uj) and a0 2 (c<j) exist except at 
isolated zeros, then we can write eq.(Bf) in terms of these cumulant generating functions as 



(u/) = 2e^cosh(^(u;)). 



(6) 



In the above ip(u) = (ipi(uj) + ^ 2 (^))/2 and ip(u>) = (ip2(ui) — ipi(uj))/2. 

It should be noted that, because of the factor a in foi^), ^2(0) = ln(a). It should further 
be noted that the zeros of <p(u), in eq.(|5|) are zeros of cosh(ip(u)) only. This is because any 
zeros of e^ cancel the poles of cosh(ip(u)). This in turn can be understood because eq.@ is 
nothing but 2(a0i</> 2 ) 5 [{0i + a0 2 }/2(a0i0 2 ) 5 ]- 

Therefore the zeros of (J){oj) may be obtained by solving 



$((jj) = ±i(l/2 + k)n = il k , 

where k — 0, 1,2, • • •. 

Now using the cumulant expansion (|2|) we can write ipiu) as 



(7) 



$(<*>) = I]7 S 



liu) 



(8) 



s=l 



where *y s = (7^ — 7<; 1 ))/2. 

If 71 = 2UIS)^ ) ^ jJ ) n 7^ ^' th e inverse function uo = ?/> _:L (■?/>) exists near the origin, and 



iu)=0 



the Taylor series about the origin gives, using the Lagrange's formula |L2 



s 

'^'=J2 b »- ! T' ( 9 ) 






where 

1 d 8 - 1 , iu 



(10) 



Notice that @ is the inverted series of (§). It should also be noticed that, since 7V s are 
real, 6 s 's are also real. This leads to theorem I. 

Theorem I. The zeros of the characteristic function 4>(u) can be expressed by 

iuj k = }_^b s — — , (11) 



5=1 *■ 



where 6 s 's Ik are given by (|10"D and (0). 

Although the dividing point x* is arbitrary and 6 s 's should not depend on x*, 7,'s have 
little meaning unless the original probability density function has two separate distributions 
over two distinct regions and x* is taken at some point between the two regions. 

For a symmetric distribution where f(x* + x) — f(x* — x), it is convenient to consider 
a random variable defined by Y = X + x* . Then the density function becomes symmetric 
about the origin, f(—y) = f(y) and the zeros of <j) y {uS) are the same as those of 4> x {uj) since 
<j> y (u) = e ix ^(j) x {u). 

It can be easily shown that the characteristic function is symmetric about the origin 
iff the density function is symmetric. Let us take for a symmetric density function f(y) 
about the origin and the corresponding characteristic functions of truncated distribution 
functions, ^(w) = 2 J^e^ f(y)dy and <j) 2 (u) = 2 / °° e*"" f(y)dy. Since ^(-w) = 2 (cu), 
we have 7, = 7^^ — 7^ = 0, for even s, and % = j^ = —jf', for odd s. 



Therefore 



(iu) 2s+1 ~ ,. u 2s+l 



m = g^(2JTI)! = Z 'S ( " 1)SWl (2JTI)!- (12) 



s=0 \"" J ' V" s=0 

By the inverse series given by eq.(||), we have 



» = ^t lYK+1 W^' (13) 

Substituting the solution for zeros (0) we have 

od t2s+1 

^B-D'W,^, (M) 

provided that series converges. Since 6 s 's are real u^s are real. This leads to the following 
corollary. 

Corollary. Zeros of characteristic function of a symmetric distribution function lie on 



the real axis, provided that the series in (11) converges 



This means that in the complex z = e~ 1 ^ '-plane, the Mellin transformation P(z), defined 
by P{z) = <f)(—i\n(z)) [|IT]] has zeros only on the unit circle. 

For some probability distributions we show next that it is possible to calculate zeros of 
the characteristic function explicitly. 

If the density function of the Gaussian distribution is f(x) = exp(— ^ )/V^cr, the 
cumulant generating function is given by ip(cu) = fi{iui) + ^a 2 (iu) 2 ||11| . Since the exponen- 



tial function of an arbitrary entire function can not have zeros in the finite region |12| and 
ip(ui) in the above is an entire function, (f)(uj) = e^^ can not have zeros in the finite region. 
This leads to theorem II. 

Theorem. II. The characteristic function of the Gaussian distribution has no zeros in 
the finite region. 



On the other hand if the density function of a double Gaussian peak is given by 

_{x- n ) 2 _(x-fi 2 ) 2 

f(x) = e 2CT i /V27TO-! + ae ^ /V2^a 2} (15) 

the cumulant generating functions for the two peaks are given by ipi(u) = fii(iu) + \a\ {ioS) 2 , 
and tp2{ u ) = hr(a) + /j, 2 (icc?) + |t7f(za;) 2 respectively. Therefore ip(uj) = ln(a)/2 + m(iu) + 
^a 2 (iu) 2 , where m = (fi 2 — A*i)/2 and a 2 = (of — crf)/2. Solving the equation (0) using the 
above ip{uS), the zeros of the characteristic function can be obtained explicitly after perform- 
ing some algebra. This leads to theorem III. 



Theorem. III. If the density function of a double Gaussian peak is given by (|15D, 
the zeros of the characteristic function are u k = Afc/|<x| + i(m — \o-I k /\ k \)/a 2 , where 
A, = ±[{(* - ^) 2 + 4¥ " (S " ^t>- 

In the asymptotic limit where a —>■ 1 and a — > 0, the zeros are given by 

u k = I k /m/(l - e k a 2 /m) + ie k (16) 

where e k = ln(a)/2m — \a 2 I\jw? . If we put a = 1 and a = in the above, then we have 
e k = 0, which makes u k s real. This is an explicit example of the corollary of theorem I. 

The above mathematical results can be readily applied to the theory of the phase tran- 
sition. As an example consider a canonical partition function Z(f3) defined by Z(j3) = 
Pf^ e~ l3E Q(E)dE, where Q(E) is the density of states at E, [3 is the inverse temperature 
l/fc&T and kj, is the Boltzmann's constant. We can identify this partition function as a 
moment-generating function A4(t), 

/oo 
e tx f(x)dx (17) 

-oo 

where the probability density function is given by 

PCX) 

f{x) = tt(x/[3 )e- x / / n(x/p )e- x dx. (18) 



In the above t = 1 — [3 / (3 , x = /3 E, and j3 is a reference inverse temperature around which 
the system fluctuates. 

Since the cumulants of f(x) are related to the derivatives of the free energy, theorem I 
implies that the zeros of the partition function can be expressed in terms of the discontinuities 
in the derivatives of the free energy provided that the first order derivative has a non- 
vanishing discontinuity If f(x) is a symmetric function, then the corollary of theorem I says 
that the zeros of the partition function lie on the unit circle in the complex e 4 -plane. Thus 
we have extended the Lee- Yang unit circle theorem to the continuum case. From the general 
principle of statistical mechanics, f(x) can be approximated by the Gaussian distribution. 
In the thermodynamic limit its mean is given by the internal energy j3 U, and the variance 
by the heat capacity as a 2 = C v /kb- Here [3 is the inverse temperature of a single phase. 
Due to theorem II, the partition function has no zeros. 

On the other hand if the system undergoes a first order transition at f3 = /3 C , then 
f(x) is characterized by a double Gaussian peak [!]|| separated by the discontinuity of the 
internal energy, or the latent heat in the thermodynamic limit. Let us further assume that 
the ratio of the weight of the two peaks is a, and there is also a discontinuity in C v , AC V . 
Let N be some integer representing an extensive thermodynamic quantity, say the number 
of particles of the system. Introducing the reduced internal energy u = j3 c U/N, the reduced 
latent heat I = Aw = u 2 — U\ and the reduced specific heat c = C v /Nki, theorem III now 
says that zeros of the partition function for small values of k may be written as, 

ln(r fe ) = U{t k ) = - ln(a)/Nl + l -{/\c/l)^ k /lf (19) 

e k = $s{t k ) = ( V0{i + \(Mi) H°)/ni 

-i(Ac/l)V fc /Q a } (20) 

Here we have used a 2 = AciV/2, m = IN/2 and & k = (1 + 2k)n/N with k = 0, ±1, ±2, • ■ ■. 
We see that the dominant finite size correction is the terms dependent on the asymmetric 
factor a. 

8 



In the thermodynamic limit where N — » oo, a dependent terms vanish, and the equation 
for the locus of zeros becomes 



r = e (A C /2 ^ (21) 

where r^ and 9k are replaced by the continuous variables r and 9. In Fig. 1, we plot u as 
a function T/T c , and zeros in the complex t— and e 4 — planes. / = 1.0 and Ac = ±0.2, 
are used in all three figures and the exact zeros given by theorem III are calculated using 
iV = 20 and a = 1. 

The angular density of zeros defined by Ng(6) = l/(9k+i — 9k) can be written for small 
values of k, as 

2^(^) = /{l + ^(Ac//) 2 ^}. (22) 

If the transition is symmetric (Ac = 0), the locus of zeros becomes the unit circle with the 
uniform density 2irg(6) = I. One should note that the equation for the locus of zeros is valid 
only near the real axis. This is because terms beyond the Gaussian approximation become 
important as the argument 9 grows, i,e for large values of k, as we have shown with an 
example in ref . || . Finially it should be remarked that the number of zeros in this example 
is infinite in the complex t-plane (Fig. 1(b)). However only a finite number of zeros closes 
the circle in the complex z = e 4 for a finite system if we consider only the first Riemann 
sheet. In fact for the symmetric case where Ac = 0, there are exactly N zeros distributed 
uniformly on the unit circle if we scale the energy of the system by /, so that / = 1.0. 

In conclusion we have shown that the scenario for the first order phase transition put 
forward by Yang and Lee || is valid in a continuum system. Becuase we have shown a 
formal relation between the discontinuities in the derivatives of the free energy and the 
distribution of zeros of the partition function, it can now be applied to any type of first 
order phase transition. It can also be extended to a multi-phase coexistence point. In this 
case, one only needs to consider the multi-dimensional complex space as we have done in 
ref. . The existence of the formal relation presented in this paper was suspected by the 

9 



original proponents of the theorem Lee and Yang themselves. In their 1952 paper 0, they 
expressed their sentiment in the concluding remark by saying, ". . . distribution (of zeros) 
should exhibit such simple regularities .... One can not escape the feeling that there is a 
very simple basis underlying the theorem, with much wider application, which still has to 
be discovered." We believe we have discovered this simple basis. 
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Korea Science Foundation through Research Grant to the Center for Theoretical Physics, 
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FIGURES 
FIG. 1. (a) it as a function of T/T c = 1/(1 — t). u\ = 0.7 and c\ = 1.0 are taken arbitrarily. 

(b) Zeros in the complex t-plane given by theorem III. For Ac > 0, the line of zeros arches toward 

the positive real axis. Only zeros between the two dashed lines appear in the first Riemann sheet 

in the complex e'-plane. (c) Zeros in the complex e'-plane. The solid lines are the loci of the zeros 

in the thermodynamic limit given by (21). For Ac > 0, the locus lies outside the unit circle which 

corresponds to the symmetric case, Ac = 0. 
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